Emulating Non-Abelian Topological Matter in Cold Atom Optical Lattices 
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Certain proposed extended Bose-Hubbard models may exhibit topologically ordered ground states 
with excitations obeying non-Abelian braid statistics. A sufRcient tuning of Hubbard parameters 
could yield excitation braiding rules allowing implementation of a universal set of topologically 
protected quantum gates. We discuss potential difficulties in realizing a model with a proposed 
non-Abelian topologically ordered ground state using optical lattices containing bosonic dipoles. 
Our direct implementation scheme does not realize the necessary anisotropic hopping, anisotropic 
interactions, and low temperatures. 
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I. INTRODUCTION 

Topological matter is operationally defined^ as a two- 
dimensional quantum many-body system with a non- 
trivial ground state degeneracy immune to weak local 
perturbations. The existence of an excitation gap sep- 
arating the ground state from the low-lying excitations 
guarantees quantum immunity of the ground state de- 
generacy against local perturbations. Topological mat- 
ter sustains, in general, two types of excitations. With 
type I (type II) topological order quasiparticles obey 
Abelian (non-Abelian) statistics. Physically braiding ex- 
citations in type I and type II topological matter mod- 
ifies the many-body wavefunction by a phase and non- 
trivial matrix, respectively. Braiding non-Abelian quasi- 
particles may enable fault-tolerant topological quantum 
computational^, reducing stringent quantum error correc- 
tion procedures required of ordinary qubit-based quan- 
tum computation. We further classify type II topolog- 
ical order based on potential application in topological 
quantum computation. Braiding a small number of ex- 
citations in type Ila topologically ordered matter does 
not yield a universal set of quantum gates necessary for 
implementing all quantum codes. Braiding in type Ila 
matter must be supplemented by unprotected quantum 
gates in order to implement an arbitrary quantum code 
thereby offering a "partial" topological immunity to weak 
local noise. Braiding excitations in type lib systems, 
in contrast, yields a universal set of quantum gates of- 
fering "full" topological immunity in implementing an 
arbitrary quantum code with braid operations. Models 
demonstrating type lib topological order incorporate ad- 
ditional complexity to accommodate a universal set of 
gates in the excitation braid structure. 

Although a great deal is known theoretically about 
topological matter in the effective field theory sense, the 
necessary conditions for the emergence of topological or- 
der and non-Abelian quasiparticle statistics in real ma- 
terials are not entirely known. A current candidate of 
type Ila topological order is the so-called ly — 5/2 frac- 
tional quantum Hall state occurring at mK tempera- 
tures in GaAs based high mobility two-dimensional elec- 



tron systems subjected to a strong external magnetic 
field. However, no direct experimental evidence exists 
establishing the non-Abelian topological nature of the 
ly = 5/2 fractional quantum Hall stated or any real 
material. Similarly, it is conjectured^ that the fragile 
and rarely observed 12/5 fractional quantum Hall state 
may be type lib with anyonic quasiparticles (the so- 
called " Fibbonacci anyons" ) suitable for universal quan- 
tum computation, but little is known about the nature of 
this very weak state. It is, however, well-accepted^-^ that 
the quasiparticle excitations of the well-known fractional 
quantum Hall states (i.e. 1/3, 2/3, 2/5, 3/7, and so on) 
are type I. 

Given the absence of any experimentally definitive 
topological system, much theoretical work has gone into 
effective theoretical lattice models which have topologi- 
cal many-body ground states. The most famous exam- 
ple of a type Ila topological lattice model in the quan- 
tum information context is Kitaev's toric codc^, a two- 
dimensional spin lattice model. Another such example, 
of interest to our work presented in this paper, is a Bose- 
Hubbard lattice model on a kagome lattice with extended 
and ring exchange interaction terms which has been pro- 
posed as a model carrying type lib topological order— i^. 
One possible advantage of the topological lattice mod- 
els is that, although these lattice models are highly con- 
trived from the solid state physics perspective (and their 
applicability to real solid state materials is completely 
unknown), it is, in principle, possible to imagine emulat- 
ing them on cold atom optical lattices similar to what 
has been already experimentally achieved in realizations 
of the on-site Bose-Hubbard model with bosonic opti- 
cal lattices^iiSiii. Motivated by the remarkable success 
of quantum analog simulation of solid state models^ 
several recent theoretical proposals have been made to 
creatoiSiii^ and manipulate^ cold atom optical lattices 
emulating Kitaev spin models^iii^. 

In this paper we theoretically identify key issues in 
optical lattice emulation of the extended Bose-Hubbard 
model originally considered as a quasi-realistic lattice 
model by Freedman et alJ^. While Refs. and d 
study this model independent of specific experiments, it 
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is natural to ask if experiments can approximate such a 
model. We analyze a "direct" optical lattice emulation 
of the topological extended Bose-Hubbard model with 
dipoles. By "direct" we mean optical lattices formed 
from interfering standing wave lasers containing bosons 
with a dipolar interaction. Our direct scheme com- 
bines recent work including optical lattice experimentsi^, 
proposals^^, and results showing Bose-Einstein conden- 
sation with dipolar atoms^. We find that our direct 
implementation scheme would be extremely difficult, if 
not impossible, even as a matter of principle because 
the combined parameter constraints on the hopping, var- 
ious interaction terms, superexchange, temperature, sin- 
gle band restriction, and chemical potential required to 
reach the proposed topological regime are, for all prac- 
tical purposes, mutually exclusive. Specifically, we find 
that: i) Our lattice setup significantly alters chemical po- 
tentials to yield prohibitively long hopping times between 
sites (See Fig. [5]) h) Band effects in our lattice setup, 
combined with the isotropic dipolar interaction, do not 
induce a sufficient anisotropy in the interaction (we re- 
quire an order of magnitude relative anisotropy among 
next nearest lattice sites but we achieve 8% at best), 
iii) The realization of many-body superexchange in har- 
monic optical lattices presents stringent constraints on 
temperature, T (See Fig. [3]). 

The plan of the paper is as follows: In section |TT] we 
review relevant aspects of the extended Bose-Hubbard 
model of Refs. and d. In section IIIII we calculate the 
Hubbard parameters for a direct implementation of an 
extended Bose-Hubbard model with dipoles confined in 
a kagome optical lattice. We attempt to modify the lat- 
tice to tune Hubbard parameters. We find that our lat- 
tice "coloring" scheme does not yield an appropriate set 
of Hubbard parameters. In section IIVI we discuss prac- 
tical issues in realizing low temperature superexchange 
with harmonic optical lattices. We find that weak (har- 
monic) site confinement places prohibitive constraints on 
the temperature. In section |V] we summarize difficulties 
we encountered in implementing the topological extended 
Bose-Hubbard model. We emphasize that our work here 
underscores difficulties with our direct implementation 
scheme using dipoles. 



II. MODEL 

We consider an optical lattice setup which approxi- 
mates a single band Bose-Hubbard model. In Refs. 7 
andd it was argued that bosons hopping in a 2D kagome 
lattice (see left panel of Fig. [1]) can realize quantum topo- 
logically ordered ground states of type lib, with Bose- 
Hubbard parameters tuned to specific values. The pro- 
posal presents the following model with an additional 
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FIG. 1: Left panel: Grey scale plot of the kagome lattice, 
defined by Eq. O in the x — y plane. Dark regions indicate 
sites. Sites encircled with solid (dashed) lines indicate red, r, 
(green, g,) sites. Sites not encircled indicate black, b, sites. 
Right panel: The same but with an additional potential (Eq. 
[6]) added to color the lattice with Is = 0.3. 



ring exchange term^i^: 



i <ij> 

+ ^^ni(ni - 1) + C/q ^ Uirij 

i (»,i)60 

+ ^ VijUiUj. (1) 
(i,j)e[>q,^0 



The fii define site dependent chemical potentials and 
is the number operator. The second term indicates bond- 
specific nearest neighbor hoppings with energy gain tij, 
where bj annihilates a boson at the site j. Uq represents 
an onsite interaction energy penalty assumed to be the 
largest energy scale in the single band limit. Uq is an 
energy penalty between particles on a given hexagon and 
is tuned to be large thereby preventing double occupancy 
of hexagons. With one particle per hexagon and uniform 
hopping one finds a degenerate manifold of boson config- 
urations that can be thought of as dimer configurations 
where each dimer lies along a line connecting the center 
of each hexagon. The degeneracy can be lifted by modify- 
ing hoppings or adding inter-hexagon interaction energy 
penalties (the last term in Eq. [T]). The latter impose 
an energy cost (to be matched with superexchange in- 
teraction energies) between next-nearest neighbors lying 
along bow-ties not within the same hexagon. A proposed 
tuning of fii, Uj, Vtj, and an additional multi-site ring 
exchange term (not discussed here) drive the system to- 
wards type lib topological order— i^. In the same context 
a simpler set of conditions were proposed for realizing 
type I ground states (the A: = 1 topological phase real- 
ized with Eqs. 9-14 in Ref. 0). Some of these conditions 
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are summarized below: 



VX/Uo 



= trb/(cC/0) 



« Uo 

t < Uo 

T < t^/Uo 

Interactions ^ A. 



(2) 



Here t^^p indicates hopping between sites colored a G 
{b,g,r} and fi e {b,g,r} with 7 e {b,g,r} the color of 
the site opposite a and (3 in the corresponding triangle, 
Fig m indicates Vij where a and /3 represent the 

colors of next-nearest neighbors i and j while 7 is the 
color of the site between them, e is a small positive num- 
ber and c and a are constants. A modification of these 
conditions, involving additional complexity, may, as pro- 
posed in Refs. J, and 8, drive the system towards type 
lib topological order. The last two conditions depend on 
temperature and the energy splitting between the lowest 
and first excited band, A. They are practical constraints 
for a direct realization of a Hubbard model in a single 
band superexchange limit with low temperatures. 



III. AN OPTICAL LATTICE OF DIPOLES 

We now consider key difficulties in directly implement- 
ing the above tight binding model using bosons con- 
fined to optical lattices. Optical lattices offer a tun- 
able environment free from defects, impurities, and lat- 
tice phonons. Implementations of Bose-Hubbard models 
have realized low temperature (T < t) superfluid and 
Mott phases^iiSiii. These Bose-Hubbard systems have 
been realized with Alkali atoms parameterized by a zero- 
range contact interaction. Effectively contributing only 
t/o in Eq. [1] Recent work seeks to extended the range of 
interaction between particles in optical lattices by pro- 
moting bosons to higher bands^iSJ^. Promoting bosons 
to higher bands uses band effects to expand the contact 
interaction to include a weak nearest neighbor interac- 
tion. In implementing Eq. [T] we need interactions to be 
tuned over several nearest and next nearest neighbors. 
Work in a different syste m^^i^^ highlights the possibility 
of confining dipoles to optical lattices. Dipoles in lattices 
generate nearest and next-nearest neighbor interaction 
terms in Hubbard models. Magnetic dipoles (e.g. ^^Cr) 
have a weak dipolar component (See, e.g., Ref. [2^ for 
estimates). For example, lattice depths yielding a hop- 
ping of Q.IEr leave a nearest neighbor interaction below 
3 X IQ'^^Er for t Pi Uq with ^^Cr. But recent propos- 
als indicate that molecules with electric dipolar moments 



may yield stronger dipolar contributions to the interac- 
tion (See, e.g., Refs. [lllli). 

To be specific we assume that bosons with a strong 
dipolar interaction can be confined to optical lattices. 
The Hamiltonian for interacting particles of mass m in 
a single particle potential defined by an optical lattice, 
VbL, is given by: 



H 



2m 



E 

k 

\Y.^int{\rk-ri\) 



(3) 



k^l 



To define a two-dimensional lattice in the x — y plane 
we assume confinement along the z-direction, Vz{z), suf- 
ficiently deep to prevent excitations out of the lowest 
confined state in the z direction. The z component of 
the wavefunction can be approximated by a Gaussian: 
(j){z) — (2/7r;^)^/'* exp {—z'^/P), where I is defined by the 
confinement frequency along the z direction. The Fourier 
transform of the dipolar interaction (excluding the con- 
tact interaction), with the dipoles oriented along the z 
direction, is then^^: 
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1 



^exp(;2fc2)Erfc(Zfcp) 



(4) 



where g is an interaction parameter, kp = {k^^ky), and 
Erfc is the complimentary error function. In real space 
this interaction is isotropic, decays as r~'^ at large dis- 
tances, and has its short range part suppressed by the 
z extent of the wavefunction. We use the above inter- 
action (including a contact interaction) to calculate the 
Hubbard parameters. We assume that the strength of 
the contact interaction and the dipolar interaction, 5, 
are independently tunable over an arbitrary range. 

The kagome lattice can be defined using six counter- 
propagating laser beams of wavevector k. By interfer- 
ing the beams at specific angles the resulting potential is 
given by. Vol = ^Vkag + V^, where^: 

Vkag(r") = X! [*^°^ (fea • r + 0a/2) 
a=l-3 

+ 2cos(fc„ -r/S-f 0„/6)]2. (5) 

With ki = fc(-I/2,-^/3/2),fe2 = fc(l/2, -V3/2), fcg = 
fc(I,0),(/!'i — —4>2 = 03 = 7'', and Iq = — i?R/2 the po- 
tential minima define an "isotropic" kagome lattice (see 
the left panel in Fig. [1]). We define = /^ma?' with 
a — Ti/k. We use the tight binding basis of localized 
states (Wannier functions) centered at sites arranged in 
a kagome lattice. To calculate the Hubbard parameters 
we approximate the Wannier functions by Gaussians in 
a variational ansatz^^. We minimize the single particle 
part of Eq. [3] with respect to four variational parame- 
ters, the Gaussian width and location, near each distinct 
site. For large lattice depths the Gaussian approxima- 
tion provides a reasonable estimate of tight binding pa- 
rameters calculated from a full band theory treatment^. 
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The large number of interfering beams used to de 
the kagome lattice (Eq. [5]) superpose to yield deep 
confine for |/o| ^ 0.5£'r. This is contrast to square o 
cal lattice geometries^ that are safely in the single-ba 
tight-binding limit for much larger lattice depths (r 
IOEr). 

For dipoles in an isotropic kagome lattice defined 
the hopping, chemical potential, and interaction 
equal distances are all uniform throughout the latt 
Note that the spatial distance between site pairs aa 
a hexagon and along a bow-tie are the same. We fi 
by direct calculation, that the requirement Uo ^ 
is therefore not satisfied by a kagome lattice of dip^ 
with a 2D spatially isotropic interaction. Furthern 
the interaction between dipoles decays as r~'^ leaving i 
rections at large distances. We ignore corrections to 
interaction beyond the next-nearest neighbor. In ph; 
with an energy gap (~ C(e^) in our example here) 
assume that the long range correction terms are mui^n 
smaller than the gap thereby preventing a phase transi- 
tion. This approximation depends on a theoretical un- 
known, the stability of the gapped topological phase with 
respect to perturbations. Emulating a topological phase 
in optical lattices would present us with an experimental 
tool to probe stability. We therefore ask if some of the 
conditions imposed on Vij (Eqs.^ can be partially met 
with the short range part of the dipolar interaction. Here 
we must match superexchange terms 0{t^ /Uq) with in- 
teraction terms, Vi^j. 

We attempt to impose anisotropy in the lattice by tun- 
ing c in Eqs.[2]with a shifted potential. Applying an ad- 
ditional set of counter-propagating beams can, at least 
in principle, be used to "color" the otherwise uniform 
kagome lattice. Recent experiments have colored a two- 
dimensional square latticei^. We have here a dual goal: 
i) Can we color hoppings in accord with Eqs. [2]? ii) Do 
band effects in the colored lattice induce substantial (or- 
der of magnitude) anisotropics in the interaction? We 
study an applied shift potential established by four ad- 
ditional beams to modify the hopping in an anisotropic 
fashion: 

Vs{r) = [cos (fe4 • (r- - -Ri)) + ^cos (fog ■ {r - Ri))^ 

+ [cos (ki ■ (r - R2)) + A' cos (fcs • (r - Ha))]^ (6) 

where = A:(l/3,0),fc5 = fc(0, \/3/6), iJi = 

a(-l/4,V3/4),i?2 = a(-l/4,5V3/4),A = 2.75, and 
A' = 1.9. The first term in the potential modifies the 
hoppings by widening specific sites in the hexagon while 
the second term maintains a balance among hoppings at 
the opposite sides of a hexagon in the kagome lattice. 
As a result of the shift the site specific potentials be- 
come modified to Vol = -^o[Mcag -t- /sl4] + Vz (see the 
right panel of Fig. [1]) . The top panel in Fig. [2] shows 
that the resulting hoppings go from uniform, at Is = 0, 
to site-dependent, at finite Is- The potential, at first 
glance, yields hoppings which approximate a tuning of 
c according to Eqs. [D Note that ^3 4 is not in accord 
with Eqs. [21 The kagome lattice in Fig. [T] appears to be 
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FIG. 2: Top panel: Hopping between site pairs around a 
hexagon as a function of the shift potential strength. Is . The 
lower left hexagon from Fig. [1] is labeled 1 — 6, where 2 is 
a green site. la approximately tunes c in Eqs. O Bottom 
panel: Chemical potential around the hexagon resulting from 
the addition of the shift potential to the kagome. 



only slightly affected but this is in fact not the case. The 
bottom panel in Fig. [2] shows that the chemical poten- 
tials are drastically modified by I4. The hopping times 
from site to site become prohibitively long with these 
large chemical potential shifts (~ En) indicating that 
weak modifications of the hoppings (30%) in this imple- 
mentation of the kagome lattice lead to unwanted large 
shifts in the chemical potential. In fact another crucial 
requirement in implementing Eq. [T] is the ability to tune 
the chemical potential locally to maintain uniform renor- 
malized chemical potentials throughout the lattice up to 
0{t^ /Uo)- We find that modifying hoppings with Eq. [B] 
leads to drastic and incompatible changes in chemical 
potentials. 

The lattice coloring scheme defined by Vs leads to 
anisotropic interactions. The colored lattice alters the 
shape of the Wannier functions in a site specific fashion 
and can therefore, in principle, induce anisotropics in the 
interaction. Eqs.[2]require that bow-tie terms, V2p^ are at 
least an order of magnitude smaller than hexagon terms, 
Uq. Since these interactions cover the same physical dis- 
tance in the kagome lattice we require anisotropics in the 
interaction among next nearest neighbors. We find, by 
direct calculation, that band effects modify the next near- 
est neighbor interaction only slightly. There is, at best, a 
8% anistropy in, e.g., the quantity (V15 — V4q)/V4q, (See 
the site definitions in the top panel of Fig. [2]) in the range 
< /s < 0.4. The variation in bow-tie terms compared 
to cross-hexagon terms was much less, < 1%, over the 
same range of Is values. 



IV. TOWARDS SUPEREXCHANGE 

We now briefly discuss practical issues in realizing low 
temperature superexchange with particles in optical lat- 
tices modeled by Eq. [31 Nearly all proposed topological 
lattice models make use of a delicate competition be- 
tween anisotropic interaction terms. Ref. |7| invokes su- 
perexchange to generate some of these interaction terms 
from an underlying Hubbard model. A variety of theo- 
retical proposals in the optical lattice setting seek to re- 
alize the equivalent of superexchange with different tech- 
niques including the use of interstitials with resonant 
interactions^^ and polar moleculesi^i^. But, convention- 
ally, superexchange in a single band Hubbard model can 
be achieved with low temperatures, T <C 0{P/U), and 
interaction strengths below the band gap, U <^ A, where 
t and U are characteristic hopping and interaction en- 
ergies, respectively. The temperature sets an absolute 
energy scale which can, in principle, be experimentally 
tuned below the superexchange limit, 0(t^ /U). In prac- 
tice, however, the realization of sufficiently low temper- 
atures in optical lattice setups remains elusive and has 
been a motivating factor in recent theoretical work (See 
Refs. HilllllSIli for example). We find rather stringent 
requirements on temperatures needed to realize conven- 
tional low temperature/single band superexchange with 
our direct implementation scheme, because, as we will 
show below, single band and low temperature require- 
ments are mutually exclusive. 

One potential solution to the temperature problem is 
to simply increase energy scales. The physical energy 
scale, Eyi, is fixed by the mass of the constituent parti- 
cles and the wavelength of the lasers defining the optical 
lattice. We exclude the possibility of adjusting the mass 
and laser wavelength. With more tunable experimen- 
tal parameters one can, however, increase t by lowering 
the lattice depth while increasing Uq with, for example, 
a Feshbach resonance. But note that the temperature 
sets a lower bound for t while A sets an upper bound 
on interaction energies. Typically, t decreases exponen- 
tially with A for large lattice depths making the last two 
requirements in Eqs. [2] difficult to realize even with ex- 
tremely low temperatures. To quantify the parameter 
window we consider one-dimensional examples below. 

The precise functional form relating t and A depends 
on details of the single particle potential defining the lat- 
tice. In solid state systems lattice sites are typically 
defined by spatially local atomic cores. The resulting 
band gaps can be quite large even for large hopping 
strengths (band widths). For comparison, consider a 
one-dimensional Mathieu lattice^^ defined by interfering 
lasers. Expanding the potential near each site yields a 
parabolic (harmonic) confinement characteristic of opti- 
cal lattices. (Note that the intensity profiles of standing 
wave laser beams will almost always yield a parabolic po- 
tential about each optical lattice site, as in Eq. O Con- 
sider the Kronig-Penney modei^S. where, by contrast, in- 
dividual sites are modeled by attractive delta function 
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FIG. 3: Hopping versus band spacing for the one-dimensional 
Mathieu (dashed) and Kronig-Penney (solid) models, Eqs. [T] 
The vertical (horizontal) dotted line indicates an indepen- 
dently tunable interaction (temperature) scale. The upper 
right quadrant of the graph indicates a range of hoppings 
and band spacings that yield an ideal limit for superexchange 
given the positions of the dotted lines. 



potentials. For well separated sites [t <C A) the band 
width in the Kronig-Penney model scales more favorably 
with band spacing: 



iKp/£^KP 



A 



3/2 



_A 



KP 



exp (^-^^/ A/Ekp 



(7) 



where £^kp = /i/Sma^p for the Kronig-Penney model is 
defined in terms of a lattice spacing aKP — a to draw 
an equivalence with En for the optical lattice Math- 
ieu problem. We derive Eqs. [7] under the assumption 
t ^ A. The square root in the exponential suggests that 
Kronig-Penney-like systems (systems with tight confine- 
ment around each lattice site) are better approximated 
by single band models over a comparatively wider pa- 
rameter range, provided equivalent energy scales. We 
show this in Fig. [3] where we plot Eqs. [71 The vertical 
and horizontal dotted lines indicate an arbitrarily tun- 
able interaction energy and temperature, respectively. 
T/E-R, = 0.005 and U/E^ = 1 are chosen as examples. 
The upper right quadrant of the graph then corresponds 
to the low temperature/single band limit ideal for real- 
izing superexchange. The Kronig-Penney hopping and 
band spacing lead to a more favorable parameter win- 
dow for superexchange. Harmonically confined sites (e.g. 
Eq. [5]) require much lower temperatures. 



V. CONCLUSION 

We have discussed potential issues entering into 
our emulation scheme of a topological extended Bose- 
Hubbard model in cold atom optical lattices of dipoles 
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finding that, even without considering the ring exchange 
terms in the proposed topological Hubbard model^, it 
will be very difficult, if not impossible, to use our di- 
rect scheme to simulate the corresponding strongly cor- 
related model of Ref. fl. We find: i) The constraint 
that the superexchange energy 0{t^ /Uo) is much larger 
than the temperature within the single-band Hubbard 
model (so that the band gap A is large compared with 
Uo and Uq) is difficult to satisfy with currently acces- 
sible temperatures in experiments, (T ~ t), using a di- 
rect emulation with harmonic optical lattices. This is 
currently a problem for most proposals making use of 
superexchange in optical lattices. A lattice with Kronig- 
Penney-like site confinement may allow a wider tempera- 
ture window. Other implementation schemesi^^i^i^, in 
conjunction with a kagome optical lattic&i^, may also be 
able to avoid the prohibitively low temperature require- 
ments, ii) Our suggested modifications to the kagome 
optical lattice, a tuning of hopping parameters with ad- 
ditional laser beams, lead to drastic and incompatible 
changes in the lattice structure itself. By tuning the 
lattice to color hoppings we find large chemical poten- 
tial shifts E^) which correspond to prohibitively long 
hopping time scales, iii) We also find that, with dipoles. 



the constraints on the anisotropic interaction are too de- 
manding for our direct implementation scheme. Specif- 
ically, cross-hexagon terms should be an order of mag- 
nitude larger than bow-tie terms, V^p ^ Uq, but end 
up comparable in our scheme, Uq- Tuning in- 

teraction anisotropy with band effects leads to only a 
small, < 8%, variation in next nearest neighbor interac- 
tion terms. A low temperature optical lattice of polar 
molecules^ may show more promise in all of the above 
categories. In spite of our somewhat disappointing con- 
clusion on the prospects of creating a topological phase 
with an extended Bose-Hubbard model using kagome op- 
tical lattices, we think that it is important to continue 
thinking about cold atom optical lattices as suitable sys- 
tems for emulating elusive type Hb topological phases by 
implementing other indirect techniques beyond the scope 
of our work. This is particularly true in view of the highly 
elusive nature of type Hb topological matter, and the fact 
that optical lattices allow for the possibility of emulat- 
ing Hamiltonians which are unrealistic (certainly in their 
pristine forms) in solid state materials. 

We thank J.I. Korsbakken, K. Park, K.B. Whaley, and 
C.W. Zhang for helpful discussions. This work is sup- 
ported by the Microsoft Q Project and ARO-DARPA. 
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